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Abstract. This paper describes a search algorithm to locate values of t where the real part of 
the Riemann zeta function, <^{cr + it), is negative for c > 1. The run-time to execute the search is 
much less than a brute-force approach and relies on certain symmetries of congruence equations 
related to the zeta function. 



1. Introduction 

The goal of this study is three-fold: 

(1) Construct an algorithm to approximate a solution to a system of congruence equations 
pertaining to the zeta function. 

(2) Use the approximations in item (1) to locate regions along the line a = 1 where the real 
part of C(f + it) may become negative. 

(3) Attempt to locate and verify instances where the real part of this function becomes negative 
along this line by a method other than trial and error. 

The value of C(s) for s = a + it in the half-plane cr > 1 is commonly given by it's Dirichlet series 

oo 

CW = E;^' Re(s)>l (1) 

n=l 

or it's Euler Product form 

c(-) = nr^W' 

p 

where p is over the set of primes. From these two expressions, little can be deduced about the range 
of Cis). However, Titchmarsh in f?", provides an analysis which gives an upper bound. We first 
take the logarithm of both sides of ^ to obtain 

oo 

logC(a + zt)--^log(l-p,r). (3) 

n=l 

For s — ao + it with (Tq > li let: 

oo 

U^{u:u^ logC(ao + *0 = - E (l - P,T"°e~'*'°s(P,.)^ 

and 

oo 

V = {v:v = <I>(ao, e,,02, ■ ■ ■) = -Y^^og {l - p^^e-^'- 

n=l 
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for independent variables {di,92,- ■ ■) with < 6'„ < 27r. Note first that U C V since tlogpn = 
0„mod 2tt. 

Using the set V, Titchmarsh placed approximate bounds on the range of the set J7 as a function 
of a for logC(s) in the half-plane a > 1. These ranges are in the form of circles with centers at 
l/21ogC(2cr) given by Q. 

[/ = |u : w e D (^1/2 log C(2a), ^ log ^) | . (4) 

Expression Q is monotonic in a, and the larger a becomes, the smaller this circular region becomes. 
For example, when a — 2, the radius of U is approximately 0.5. This means the values of log C(o'+ii) 
for cr > 2 are all (approximately) contained in a circular region of radius .5 and center at (0.04,0). 
For all values of u > 2, the absolute value of logC(s) is smaller than 0.504 < tt/2. Consider ((s) = 
—a+i7 for positive a somewhere in the half-plane cr > 1. Then log(— 0+47) = hi\—a+i^\+i{6+2mT) 
with \9 + 2n7r| > 7r/2. But tt/2 > 1/2 and so ReC(c -f- it) is never negative for a > 2. In order for 
the real part of (^{a + it) to be negative we must have 

or cr < 1.197. Once a becomes smaller than this value, the Titchmarsh circles extend beyond the 
tt/2 line and so ReC(s) can and does become negative. For the purpose of this study, we only 
investigate C(l + ii) since by continuity, if the function has real part negative on the line cr = 1, 
then it is also negative for some neighborhood 1 < cr < 1 + e. 



Consider the expression: 



2. Initial Considerations 

N 

-^log(l~p-'^e-'*) (5) 

n=l 

for a constant cr. What value of t causes the fewest terms of this sum to reach an imaginary part 
outside the interval (~7r/2, tt/2)7 It's easy to see that one such value ist ^ tt/2 since in that case, the 

terms in ^ become log ^1 + and thus give rise to the largest argument. Table jlj summarizes 

how this build-up of arguments proceeds. By the time N — 14, the imaginary component has 
exceeded |7r/2|. Now consider the partial sum: 

N 

-5]log(l-p-ie-^"°s(P„) 

n=l 

When p = 2 and t — ^^^^^^^ for example, the imaginary part of this term is at its maximum value 
or arctan(l/2) w 0.463. And in the case of a value of t which simultaneously reduces t\og{2) 
and t\og{3) mod 2 to a value close to 7r/2, the net imaginary part of these two terms would be 
approximately 0.79. This is the case for example, when t — ]p^p|- And when t is larger still, 

additional " confluences" of this sort can push the first n terms of the sum even closer to a negative 
imaginary part. 

These considerations do not guarantee the infinite sum will have an imaginary component outside 
the range {—tt/2, tt/2) but does give us some indication how the real part of zeta might grow negative 
in this part of the half-plane: a value of t which when reduced modulo 27r for the first n terms, 
builds up to an imaginary part greater than 7r/2 that is simultaneously not offset by the remaining 
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terms in the sum. The algorithm described in this paper searches for values of t which maximizes 
this build-up of initial argument for the first n set of primes, then searches this set of points for 
values that give rise to a negative real part of ({1 + it). 

3. Description of the Algorithm 

All algorithms for this study were written in Mathematica 7.0. The code locates values of t 
which approximately satisfy (|6| for a set of primes within a tolerance of S. One might initially 
guess that this would reduce to nothing more than a brute-force search for values meeting the 
tolerances. However, the modular space of ^ contains a marvelous symmetry that allows us to 
locate approximate solutions to ([6| in much less time than a brute-force search. The solutions to 
^ in this paper are called "confluence points" and denoted by C„(6', 5) representing a solution to 
the first n equations of ([6| for particular values of 9 within ±S. In this paper, sometimes when 
referring to these points, only C„ or "confluence order" is used or if 9 is implied, only Cn{S) is 
stated. 

'ilog(2) =£i(f,(5) mod (27r) 

Mn=i : (6) 

t\og{pn) ^D{^,5) mod(27r) 

A value of t which satisfies the first n equations in ^ will push the imaginary term of the log sum 
into a region outside the interval (— 7r/2, 7r/2). If the net effect of the remaining terms of the sum 
do not appreciably build-up in the opposite direction, then this value of t will cause Re ^ to become 
negative. The plan then is to locate a large number of C„ points and check the value of C(l -1- it) 
in a neighborhood of each to determine if its real part is negative. We search along the line a = 1 
since the greater a is, the smaller the argument becomes for each term in the series and the less 
likely the initial terms will have any affect on the final sum. 



Table 1. Table (defauh) 



N 


Inilog(l + ipn) 


1 


0.46365 


2 


0.78540 


3 


0.98279 


4 


1.12469 


5 


1.21535 


6 


1.29212 


7 


1.35088 


8 


1.40346 


9 


1.44691 


10 


1.48138 


11 


1.51363 


12 


1.54065 


13 


1.56503 


14 


1.58829 


15 


1.60956 
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Figure 1. Modulus space M, 



4. Mapping the modular space of M„ 

The algorithm described below is designed around a "base confluence set". We start by deter- 
mining the confluence points for two primes. The following description explains a base set in terms 
of the first two primes. The algorithm is conceptually the same for at least the pair (5, 7) and likely 
so for other pairs. 

Consider the modular equations 

7712 (i) = mod (tlog(2),27r) - 7r/2 
7713 (t) = mod (Hog(3),27r) - 7r/2, 

the plots of which are show in Figure [T] with 77i2(i) in blue and ra^{t) in red. Where the two lines 
cross, we have a solution to the system: 

<log(2) = .T mod (27r) - 7r/2 

tlog(3) = X mod (27r) - 7r/2 

for some — 7r/2 < x < 37r/2. One immediately notices a symmetry in Figure [T] the confluence 
points, shown as black dots, line up. If we connect the dots with diagonal black lines running from 
the upper left corner to the lower right corner and expand the range of t, we obtain the plot in 
Figure [2] (to avoid clutter, the modular 7n„ lines have been removed). Equations for the black lines 
are given by 

L„(t) = (31og(2) - 21og(3))t + 27i7r-7r/2, 71 = 1,2,3, ••• . 

The lines Ln(t) serve as a first approximation for locating the next order of confiuence points: we 
simply calculate where the black lines cross the t-axis, check if a base confluence point is near-by, 
then check for a C3 match within a specifled tolerance. However, the L„ lines are too-closely spaced 
to offer an efficient algorithm for locating a large number of C3 points quickly. We can improve the 
algorithm by considering a second set of symmetries in Figure [2] which is shown as the green lines 
represented by the equations: 

Gnit) = 1/3(8 log(2) - 51og(3))i- 27i7r/3-H37r/2, 71 = 1,2,3, - •• 

In order to further reduce the number of points that need to be checked, we can again form a 
system of diagonal lines connecting the base confiuence points. These are shown as the red lines in 
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Figure 2. Diagram of L„ and G 
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Figure |3] where this time, both the mod Hnes m„ and black hnes L„ were removed for clarity. The 
red lines are given by the equations: 

R{n, m, t) = log (4/3(2/3)22/31^ t + 7r/62(16n + 4m - 47) 

for m — 1, 2, 3, 4 and n = 0, 1, 2, • • • . And although we could continue this process to obtain even 
further improvements, we stop at the R{n,m,t) lines for this study. Note in Figure |3j the R lines 
cross the i-axis with a very small slope with the two nearest C2 confluence points along each red 
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line being very close to the f-axis. Where each R lines crosses the t-axis, the nearest-neighbor set 
denoted by ^(2,3), satisfy 



B r^/9A^ /Hog(2)^i^(f,5) mod (2^) 
^'^'^'^"/''^^ = \aog(3)^I?(f,^) mod (2.) 



for small S. An example nearest-neighbor set for one R line is show in the figure as the two red 
dots. It is not difficult to determine the zeros of the lines R{m, n, t) and then subsequently,the 
two nearest neighbors. This allows for an efficient search of the modular space for the points 
C2(7r/2, 5) requiring only a few computations for every 2000 increment in t. Each confluence point 
C2, along the R lines is at a distance of A^(2,3) = iog(3)-Tog(2) ~ "^^^ distance between 

each successive zero of the R lines is 53iog(3)^s4 iog(2) ~ 3009 and the slope of each R line is 
thr = log(4/3(2/3)22/3i) « -6.7x10"^ Thus, knowing the zero point of each R{m, n) and A^(2, 3), 
we can determine with relatively few calculations, the nearest C2 confluence points closest to the 
t-axis and then check only these for higher-ordered confluence points. But it turns out that we can 
do much better than this because the pattern of confluence points in one range of values of t goes a 
long way in determining the pattern in a following range of values. Some drift is of course expected 
but in actual runs of the software, this turned out to be surprising little over a very large interval. 
We can exploit this property by using a set of lower-ordered confluence points as a map for the 
next higher set of confliience points thereby greatly reducing the time to locate these points. 

It's easy to see that the maximum height reached by either the previous or next C2 point along 
each Rn line is mn iog(3^^^og(2) ~ 0.032. Therefore, if the tolerance 5 is set greater than 0.032, than 
both the previous and next point will be below this value and can be chosen without checking the 
residue mod 27r. 
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The following Mathematica code implements this approach to locate the base- confluence points 

-8(2,3)- 

For [n=l ,n<nMax,n++, 
For[m=l,m<4,m++, 

(-31 + 4m + 16n)7r 



zeropt = 
startpt = 



2(531og(3)-841og(2))' 

2(31 - 7m + 3n)7r^ 
log(3) - log(2) ' 



62 

intpt = IntegerPart[— (zeropt — startpt)]; 

k 

( 627rinpt \ 
startpt H 1 ; 

/ 627r 
nextpt = I prevpt + 



k 

If [theDegree<maxDif f erence , 

thePoiiitList23=Append[thePointList23 , prevpt] ; 
thePoiiitList23=Append[thePointList23, nextpt] ; 

> 

testprev=Mod [prevpt Log [2] ,27r]-7r/2; 
testnext=Mod [nextpt Log [2] ,27r]-7r/2; 
If [0 < testprev < theDegree, 
thePointList2=Append[thePointList2, prevpt] ; 

]; 

lf[0 < testnext < theDegree, 
thePointList2=Append [thePointList2 .nextpt] ; 

]; 

]; 

]; 
]; 

Prom this starting set of points, wc go on to flnd the set of C3 points by calculating the gaps 
between successive C2 points and then use those gaps to predict with good success, the location of 
the next C3. In particular, if we use the above code to locate the first 10, 000 C2(7r/2, 0.01) points, 
we find only three possible gap sizes: 

r 7787r 73607r 81387r 1 
32, 53} - I j^^^, j^^^, j^^^l 

We can then use this table to search for C3 point by simply checking three possible locations. For 
example, if a C2 point was located at t2, then we need only check the three points ^2 + 5i, ^2 + 52) 
t2 + g3. The following code implements the next stage of the search by locating the C3 points 
using the previously-calculated C2 gap table. In the code, we have a While loop which checks 
for the C2 points against the corresponding gap table. As the confluences grow in order, we add 
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conditionals for additional primes. The following code is only the main routine. In actual practice, 

pre-processing of the previous confluence data is necessary. The variable mintval is then initially 
set to the last C„ point (if any) found in the pre-processed data. After each acquisition, the data 
must then be post-processed into a sorted gap table for the next confluence search. 

thePointList3 = {}; 
mintval=last in23; 
n=l; 

For inypoints=l,inypoints < iiiaxN,mypoints++ 
n=l 

testpt = mintval -|- gaptable23[[n]]; 
While [ 

Not[0< Mod [test Log[2] ,27r] -7r/2 <theDegree] , 
n++; 

If [n>gaptable231en, 
Print ["Gap table 23 exhausted"]; 
Abort [] ; 

testpt = mintval -|- gaptable23 [ [n] ] ; 

]; 

]; 

mintval = testpt; 

If [Mod [mintval Log [5] , 27r] — 7r/2<theDegree , 
thePointList3 = Append [thePointList3, mintval] ; 

]; 
]; 

With this code and some additional pre-processing and post-processing code, we are able to find a 
large number of C^{'n /2^5) points quickly. For example, we can check about 50,000 C2 points in 
about 40 seconds with the machine this code was tested on. We next calculate the gap table for a 
sufficient number of C3 points. In actual practice, this was usually around 1500 points. In one run 
of the C3 code, we obtained the following gap table: 



_ r 9126077 3600147r 4512747r 9219367r 10131967r 13732107r 14644707r 1 
gaptable3 - \^^^^y fo^(372) ' fo^M ' log(3/2) ' log(3/2) ' log(3/2) / 

Notice the gap table has grown in size from the first table. This is quite expected considering we 

are working with residues mod 2it. The next step of course is to use gaptable3 to locate the 
C4(7r/2, 5) points. We then add another conditional to the While loop to check for a C3 confluence 
point. We then create a gap table of those points and continue the search and if necessary, adjust 
the tolerance. The routine is coded to terminate a run if after checking all values suggested by the 
gap table for a particular point, the routine fails to find the desired confiuence point at the specified 
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tolerance. In this case, a small change to the tolerance was found to correct for the drift at least in 
the range studied. 

5. Acquisition of Data 

All code was run on a 2.2 GHz machine with one Gb RAM. In order to study this algorithm, 
three runs were performed. The first two were run with a base-confluence set i?(2.3)- The details 
of constructing this set were described in the previous sections. The third set was done with a 
base-confluence set i3(5 7), and it's construction is virtually identical to the S(2,3) set except the 
equations for L„, Gn, and i?„ are different reflecting the use of different primes in their construction. 
In all cases, 9 = ^. 

5.1. Run 1. Using base-confluence ^(2,3): and S set to 0.1, the algorithm located approximately 
1000 points for each confluence up to Cio in approximately 90 minutes. The variable nMax was 
set to 50,000. This means that 50,000 C„_i terms were checked in order to obtain the C„ set 
within the currently set tolerance. In all cases, the tolerance was "globally" set meaning that each 
modular calculation was checked against the same value of 6. Table [T] summarizes the result of 
this run. "First" and "Last" are the first and last points of the set, "Gap" is the average distance 
between C„ points as determined by this algorithm, and k — iog(3)!!^iog(2) 1 that is, the first C3 point 
is located at ^^(fP^. 

Table 2. Run 1 



Type 


Points 


First 


Last 


Gap 


^(2,3)(-l) 


7992 


432fc 


1551976fc 


10^ 




923 


8274fc 


9342654fc 


10^ 


C4(.l) 


817 


171406fc 


450025676fc 


10^ 


C5(.2) 


1601 


6538724fc 


11407973800fc 


108 


C6(.2) 


1625 


180883380fc 


101681208846fc 


109 


Cr(.2) 


1659 


1524427526fc 


556034603992fc 


109 


C8(.2) 


1689 


34995876276fc 


3927590605616fc 


lOi" 


C9(.3) 


2438 


232692763660fc 


20585726758472fc 


10" 


Cio(.4) 


3197 


5206151247198fc 


53443637744022fc 


10" 



Note the tolerance was relaxed at C5 and Cg due to gap table exhaustion. Figure |4] represents a 
"confluence portrait" for one Cg value obtained in this run and serves to graphically illustrate the 
solutions for Afg of ([6]). Each line in the figure represents one equation of ([6| plotted in the range 
(< — a, ^ + a) for some a. Notice how all the equations confluence at this Cg point hence the name 
given to these points. 

Consider the six terms of the partial sum of ([5| which for simplicity we take a — \ since we are 
only interested in an approximation: 



6 

E 

n=l 



log 1 



,-lt6 log(p) 



V 



-0.126968 - 1.32215i, 



and note the imaginary part is already close to |7r/2| and even although this gives no indication of 
how the remaining terms will contribute to the sum, it illustrate how the confluences are at least 
affecting the first few terms of the log sum. 
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Figure 4. C^i-K /2,Q.2) at U = 20 430 730 768/c 




Figure 5. Cio(7r/2, 0.08) at t = 34494 360 555 864 694fc 

5.2. Run 2. The second set of data was initially acquired at (5 = 0.05 (except the base-set) and with 
a sampling size of 200,000. The gap tables successfully located on average 1600 confluence points 
until Cg was reached. Then 5 was relaxed to 0.07 and the Cg set successfully acquired. Alternatively, 
we could have re-acquired a larger gap table for Cg. The acquisition of Cio exhausted the Cg gap 
table. This set of confluence points was re-acquired dX 6 = 0.08 until after the acquisition of 1273 
terms, the gap table was again exhausted. Total execution time was approximately five hours. 
Figure [5] shows a typical Cio confluence portrait at this value of 5. Note the equations at the origin 
confluence at a smaller tolerance compared to those in Figure |4] 

5.3. Run 3. Since one objective of this study was to devise an efficient algorithm for numerically 
solving the system of congruences defined in one would like to know how well it worked for 
extremely large values of t. This third run attempted to determine this by running the algorithm 
at an initial tolerance of (5 = 0.01 and sampling size of 500, 000. As Table |4] suggests, the algorithm 
was stable enough to acquire 1765 C14 points with an average gap interval of 10^^ at a tolerance of 
0.05. This allowed the acquisition of 121 C15 before exhausting the C14 gap table which by then 
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Type 


Points 


First 


Last 


Gap 


5(2,3)(-01) 


1239 


3186fc 


1551556fc 


10-* 


C3(.05) 


1599 


118922fc 


24222G222fc 


10^ 


C4(.05) 


1602 


66656806fc 


19751098706fc 


108 


C5(.05) 


1642 


1095953254fc 


1684031570650fc 


IQio 


Ceim) 


1624 


10971337166fc 


46958108569642fc 


10" 


Crim) 


1563 


5837564637802fc 


963484780912880fc 


1012 




1608 


2354G3791542126fc 


41141345684941152fc 


1014 




2199 


1409137033297936fc 


838179310625884816fc 


1015 


Cio(.08) 


1273 


3449430555864694fc 


5206649605555820266fc 


1016 



had grown to about 2000 entries. This data took approximately 24 hours to acquire. The final 121 
Ci5 points were checked for a minimum Cig confluence which was found at a tolerance of 0.063 and 
is shown in Figure [61 In Table lil fc = wfyil^io^fs-) • 



Table 4. Run 3 



Type 


Points 


First 


Last 


Gap 


5(2,3)(-01) 


835 


1188fc 


105918fc 


10^ 


C3(.01) 


239 


415346fc 


188449608fc 


10^ 


C4(.01) 


1760 


951398fc 


341342656fc 


10^ 


C5(.01) 


785 


945400159026fc 


248625146191138fc 


1012 


C6(.01) 


782 


358872295235106fc 


53873791685357248fc 


1015 


C7(.01) 


2449 


25465934548192130fc 


24901118319251218470fc 


10" 


C8(.01) 


1607 


4228669878993511528fc 


3026698278169709217444fc 


1019 


C9(.01) 


1581 


443049036634051704488fc 


106247940250479965840452fc 


1020 


Cio(.02) 


1608 


15784886571659169335556fc 


2201816585968897375148202fc 


1022 


Cii(.03) 


2418 


325647381902527914552536fc 


37724324922892546494833526fc 


1023 


Ci2(.03) 


1800 


4833028477320576769793928fc 


502072051726121146180272476fc 


1024 


Ci3(.04) 


3924 


60624317376187697337278868fc 


5796475948831968271200664166fc 


1025 


Ci4(.04) 


1765 


510082174019371223753727784fc 


19898784954288834658068768846fc 


1026 


Ci5(.05) 


121 


6450267346766950732885822366fc 


26894796143731914838108040998fc 


1027 



6. Processing the Data Sets 

For the first run, we tested each entry in the confluence tables up to Cg by calculating C(l+i(i±l)) 
at 20 equally-spaced points in the interval and obtained a total of ten regions along the line Re(s) — 1 
for which the real part of zeta became negative. Sample points in the indicated C„ region are given 
111 Tabled The value of zeta in the regions was determined using the Mathematica " N" command 
using arbitrary precision arithmetic with the flnal answer reported to thirty digits of accuracy to 
the right of the decimal place using the following command for each of the ten points in Table [5} 

N[zeta[l+it] ,{oo,30}] 
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Figure 6. Ci6(7r/2, 0.063) at t = 12253527959225463513391519458A; 

The smallest value of t detected in the data sets described above in which ReC(l + it) < was 
found in the vicinity of 3.4 x 10^ as shown in Table [5] (recall the factor k). All the values in Table 
[5] had negative imaginary part due to the search criteria selecting 7r/2 residues. No points in the 
first four confluence tables of the remaining runs were found to give rise to a negative real part of 
C at least using the method described above. The remaining confluence tables were not inspected 
at this time because of the impracticality of actually computing the zeta function (on the machine 
used to analyze the data) at the relatively large values of t in these tables. 



Table 5. Points where ReC(l + it) < 



ID 


Type 


t 


C{i + tt) 


1 


C3 


4 378 640fc 


2/5 


-0.009- 


1.22i 


2 




415 782 314fc - 


2/5 


-0.024- 


1.23i 


3 




20 430 730 768fc- 


1/10 


-0.015- 


l.OSi 


4 




25 705 015 862fc- 


1/5 


-0.002- 


1.20i 


5 


^6 


47 668 373108fc- 


1/10 


-0.027- 


0.96i 


6 




53 761507682fc 




-0.009- 


l.OOi 


7 




62 484 882 686fc- 


1/5 


-0.013- 


i.m 


8 




65 421460 042fc- 


1/10 


-0.003- 


0.91i 


9 




97190 286104fc- 


1/10 


-0.018- 


l.OOi 


10 




99154 858182fc- 


3/10 


-0.026 - 


1.25i 



7. Conclusions 

One can easily locate a value of t where ReC(l + it) < by trial an error by simply sampling 
the t-axis at intervals of 1/10 starting at t = 0. Doing so, one first finds a negative real part at 
t = 682112.9. However, this study attempted to predict where the real part of the function would 
turn negative without a trial and error approach. A surprising result to come from this study is 
the relative stability of the method. One would have guessed that such a method would quickly 
degrade due to chaotic drift in the modular space rendering it virtually impossible to predict any 
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future values let alone values up to 10^^. It would be interesting in a further study to try and 
explain what is the source of this stability. 
Additionally: 

(1) What are the precise shapes of these negative contours where KeC{o' + it) < and how do 
they change as we go up the i-axis? 

(2) This study focused on two base-confluence algorithms: (2, 3) for the first two trials, and 
(5, 7) for the third trial. It is not known how the algorithms are affected by different base 
prime pairs. 

(3) This study demonstrated some structure in the modular space of (|6|. Is there a similar 
periodicity in the domains where Re C(cr -f it) dips below the a — t plane? Can one predict 
where such domains might be located analogous to how the various confluence points were 
"predicted" using the gap tables? 

(4) This study did not attempt to investigate in detail the drift of values in the modular space. 
When a gap table became exhausted, the tolerance was simply increased or a greater number 
of Cn-i values were located and a new larger gap table created. It would be interesting to 
better understand this drift. 
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